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Light propagation in a medium made of densely packed dielectric spheres is investigated by using 
a rigorous diffraction theory. It is shown that a substantial suppression of the local density of 
states occurs in spectral domains where the single constituents exhibit Mie resonances. The local 
density of states decreases exponentially at the pertinent frequencies with a linearly increasing spatial 
extension of the aggregated spheres. It is shown that a self-sustaining random arrangement of core- 
shell spheres shows the same fundamental characteristics. Such approach offers a path towards easy 
to fabricate photonic materials with omnidirectional gaps that may find use in various applications. 

In the past several years light propagation in random media composed of discrete scattering objects has attracted a 
considerable deal of interest. It derives its fascination not only from numerous potential applications such as, e.g., the 
random laser [U[2], but also because it challenges our fundamental understanding of light propagation. Most notably, 
the effect of light localization attracted considerable attention [3 . Light localization can be understood as a synonym 
for the occurrence of a complete photonic band gap where spectral components of light within the photonic band 
gap must not propagate within the medium because Maxwell's equations provide only evanescent wave solutions. By 
placing a source that emits light at a frequency within the gap inside such material, the light transport off the the 
source is suppressed. It remains localized. Light localization in one-dimensional structures may be easily achieved 
by, e.g., introducing a defect into a sequence of dielectric layers arranged to form a Bragg stack. In two-dimensional 
media made of infinitely extended cylinders of high permittivity this localization was feasible by exploiting scattering 
resonances of sufficient strength [H [5] . As a signature of localization one may use the exponential decay of the local 
density of states (LDOS) at a point inside the medium with increasing size of the system. A significant suppression 
of the LDOS in a finite system may be interpreted as a sign of the appearance of a complete band gap, although it 
does not constitute a rigorous proof. Because scattering resonances also affect the formation of a gap for a periodic 
arrangement of the cylinders, a complementary computation of the photonic band structure proved that the band 
gaps occurs indeed in the same spectral domain as the reduced LDOS [5 . 

In three-dimensional systems, made of spheres [6 , this possibility of light localization is controversially discussed as 
in experiments absorption may potentially shade localization effects . To investigate theoretically the possibility of 
a suppression of light propagation, one may take advantage of an analytical approach [8 . But as the structure involves 
a high index contrast and a size of the scatterers that is comparable to the wavelength, usually the problem requires 
a rigorous solution of Maxwell's equations [9l[T0]. Using a super cell approach rigorous computations of the photonic 
density of states can be performed. For an amorphous photonic material consisting of a continuous-random-network a 
complete photonic band gap has been identified recently [11 . This work has shown that appropriate short-range order 
might also permit for light localization, in addition to the well-known suppression of light propagation in periodic 
media due to Bragg resonances [12]. 

However, most experimental investigations of light propagation in 3D random media [13 employed sufficiently 
dense packed high permittivity spheres. This is mainly due to their easy fabrication procedure that relies on colloidal 
chemistry. The excitation of Mie resonances in these spheres is usually regarded as the primary reason for suppressing 
propagation [HI [15]. To theoretically verify the effects associated with light localization, various efforts were under- 
taken in the past. In Ref. [16] the density of states was computed for small clusters of spheres using a generalized 
Mie theory. In Ref. [171 scattering effects depending on the cluster size were studied for low index contrast using the 
superposition T-matrix method. A large scale finite-difference time-domain (FDTD) method was also used to trace 
light propagation in ensembles of spheres with a strongly varying size [T8] . 

In this work we focus on investigating the LDOS of large clusters made of randomly arranged identical high 
permittivity spheres. For this purpose we use the FDTD method [19^. It is shown that the LDOS decays exponentially 
with increasing size of the cluster in certain frequency intervals indicating a significant reduction of light transport 
velocity. The spectral positions of these domains are correlated to Mie resonances. They occur independently of the 
arrangement and are tunable by modifying the permittivity of the spheres. However, below a certain contrast the 
frequency intervals of the suppressed LDOS are hardly to recognize. We also investigate the LDOS of core-shell spheres 
which were arranged by using the sphere dropping method. Consequently they form a self-sustaining ensemble of 
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spheres. The approach ahows to maintain the main optical characteristics if the permittivity of the sheh is sufficiently 
low. This work opens a path towards an easy to fabricate photonic material that allows for a suppression of light 
propagation regardless of its direction and polarization. Such a material can be applied as, e.g., isotropic filters for 
solar cells. 

In a first investigation we applied our procedure to a system in where the spheres are arranged on a diamond lattice 
[20] . The purpose is to verify our computational approach as well as to investigate the appearance of the LDOS in the 
presence of a complete photonic band gap. The results are compared with photonic band structure computations of 
the infinite lattice, proving the presence of a complete band gap. We have chosen the spheres to have a permittivity 
of e — 11.56 and a radius of = 0.25a, with a being the lattice constant. The finite diamond lattice was defined by 
a 7a X 7a X 7a cube. The spatial permittivity distribution was discretized in the FDTD method with a resolution 
of 0.025a. To compute the LDOS, we located a point source at the center of the structure and computed the time 
dependent transmitted power through a sphere that includes the entire structure. The point source was excited 
by a temporal Gaussian pulse. Fourier transformation of the time dependent signal yields the frequency dependent 
transmission. Upon a normalization with the power spectrum of the source, we obtain the transmission. In a lossless 
medium this transmission corresponds to the LDOS associated with the spatial point where the source was placed. 
On a logarithmic scale the spectrum of the LDOS as a function the normalized computational time is shown in Fig. [l] 
(a). 




FIG. 1: (color online) (a) Frequency dependent LDOS on a logarithmic scale for a finite system of dielectric spheres arranged 
on a diamond lattice as a function of the computational time, (b) Comparison of the band structure for the respective infinite 
periodic lattice with the LDOS of the system. The significant reduction of the LDOS clearly corresponds to the spectral region 
of the complete band gap. 



After a sufficient computational time of ~ 100 Aa~^ all transient effects are faded away. The LDOS converges 
to a stationary distribution. The logarithmic scale clearly reveals a strong suppression of the LDOS in a frequency 
interval between 0.4 and 0.44aA~^. The main features of the LDOS are independent on the exact position of the 
point source within a unit cell. Comparing the LDOS with the band structure (Fig. [l] (b)) clearly reveals that the 
spectral region of a suppressed LDOS coincides with the band gap. For a spatial extension of the structure of 7a in 
all dimensions, the LDOS is decreased by approximately two orders of magnitude when compared to spectral regions 
for which propagation of light is allowed. 

After having identified how band gaps and reduced LDOS are related, we proceed in analyzing random arrangements 
of the same spheres where we keep the filling fraction of the spheres constant at 41%. Although a ceases to have 
a meaning of being the period we use it further as the length scale. Exemplarily, Fig. 2] (a) shows a particular 
arrangement of the spheres. The arrangement was achieved by generating at first randomly a coordinate for a new 
sphere. Subsequently it was verified if a sphere is already present in the volume to be occupied by the new sphere. 
If this is excluded, the position is accepted, otherwise it is rejected. The procedure is repeated until the predefined 
particle density is met. The exclusion of penetrating spheres was necessary to avoid that resonances associated with 
the spherical shape of the particles (Mie resonances) were appreciably altered. Moreover, as a result of preliminary 
numerical experiments, nearly touching spheres had to be excluded. This was ensured by defining the forbidden volume 
slightly larger. Nearly touching spheres are sufficiently strong coupled to lift the degeneracy of the spectral position 
of a Mie resonance, causing an inhomogeneous line broadening and a damping of the resonance of the ensemble. As a 




consequence, the LDOS would not have been sufficiently suppressed. Hence placing a new sphere is inhibited within 
a radius of 0.3a of an existing sphere. This seemingly limitation will be subsequently lifted by covering the spheres in 
the simulation with a core material having a low dielectric constant. Such arrangement permits for a dense packing 
of the spheres but maintaining the characteristics of the LDOS. 

The LDOS of a medium composed of randomly arranged spheres as a function of the spatial extension of the 
domain that hosts the spheres is shown in Fig. [2] (b). At low frequencies the LDOS remains constant, indicating 
that light flows off the source. The source is always located in the center of the spatial domain of interest. On the 
contrary, at slightly elevated frequencies two bands exist were the LDOS shows the exponentially decaying strength 
with increasing linear size of the system. It suggests that the transport velocity of light is strongly reduced. The 
two bands are centered at ^ 0.59 aA~^ and at ^ 0.81 aA~^, respectively. By using different random arrangements 
slight modifications in the LDOS were encountered, though the principal features of a suppressed LDOS in the same 
spectral domains was always observed. 

Both peculiar frequencies can be correlated to Mie resonances of the isolated dielectric sphere [2J^ . For this purpose 
the scattering cross section is shown in Fig. |3](a). Two sharp resonances can be seen at 0.56aA~^ and at 0.81aA~^. 
From the Mie coefficients it can be deduced that they correspond to resonances of the lowest eigenmodes of the electric 
(co) and the magnetic (do) potential inside the sphere, respectively. The reason why the suppression of light transport 
occurs at slightly larger frequencies compared to the Mie resonances lies potentially in the requirement that the field 
scattered by the sphere oscillates tt out of phase with respect to the incident field. This will effectively suppress 
the propagation of light in the forward direction since scattered and illuminating fields annihilate due to destructive 
interference provided that the scattering strength of either sphere is sufficiently strong. This required phase difference 
is encountered, as in all resonance phenomena, slightly above the resonance frequencies. 

Modifying the permittivity of the sphere will have various implications. Foremost, the spectral position of the 
resonances shifts towards larger frequencies, the resonance gets weaker and broader, and overall, the scattering 
strength is reduced. These modifications of the LDOS for randomly arranged spheres are shown in Fig.|3](b). The 
spatial extension was chosen to be (7a)^. The spectral regions of the reduced LDOS shifts to higher frequencies being 
in perfect agreement with prediction of Mie's theory. Furthermore, the anticipated smaller reduction of the LDOS 
for a smaller permittivity is evident. If the permittivity of the spheres is smaller than e 6, defining a domain of 
suppressed LDOS is pointless. 

Because the structures investigated up to now cannot be fabricated, we now investigate core-shell particles. The 
shell is made of a material with a sufficiently low permittivity. It should not affect the strength of the scattering 
resonance of the core. Furthermore, the shell should have an appropriate thickness to allow for a sufficient separation 
of the core spheres. As mentioned we have chosen a core-shell radius of 0.3a with a shell permittivity of £ = 2.25. The 
core features were kept invariant {R = 0.25a, e = 11.56). The self sustaining ensemble, reasonable for fabrication. 
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(^) Frequency in ax'^ (b) Frequency in ax 



FIG. 3: (color online) (a) Scattering cross section calculated using Mie theory for a single sphere of radius R = 0.25a and 
permittivity s = 11.56. The sphere is embedded in air and illuminated by a plane wave, (b) LDOS for randomly arranged 
spheres as a function of the sphere's permittivity. 



was generated with the sphere dropping method. In this method a sphere is dropped into a fictitious box. Gravity 
forces the sphere to sink down (along the 2:-axis) to the bottom of the box. If the dropped sphere touches a another 
sphere a translation in the x — y-plane of the dropped sphere is performed. The translation is chosen such that the 
distance between the dropped sphere and the other spheres is maximized. It allows to sink the dropped sphere further 
down until it touches again another sphere. The procedure is repeated until the dropped sphere cannot sink further, 
although it is incrementally replaced in the x — ?/-plane in an arbitrary direction. Once this stage in the procedure is 
reached, the sphere is stably trapped by its surrounding. 

The procedure allows to generate self-sustaining ensembles of randomly arranged spheres with structural integrity. 
To force the randomness we assumed hard boundaries. The final spatial extension of the structure is the same as 
before and amounts to (7a)^. On passing we note that in small spatial domains a certain periodicity is encountered 
as the spheres tend to form an opal structure. But as no energy minimum for the global arrangement was forced, the 
structure will retain its random character in general. A distribution of spheres generated with this method is shown 
in Fig. [4] (a). The filling fraction amounts to 54%. 

The corresponding LDOS for such a core-shell structure is shown in Fig. [4] (b). For the sake of comparison we 
show also the LDOS for the structure without shells. In this case the larger density of spheres results in a stronger 
suppression of the LDOS when compared to the scenario analyzed before. However, the spectral position of the 
suppression dips coincide. The optically weak shell causes primarily a slightly increased LDOS in the spectral domain 
of interest, though it remains some orders of magnitude smaller compared to spectral domains where light propagation 
is allowed. A slight shift of the spectral domains of reduced LDOS appears for the second resonance frequency of 
0.805aA~^. Thus we conclude that the shell does not significantly affect the properties of the photonic material, 
though a minor degradation can be observed. On the other hand, the shell is sufficiently thick to guarantee the 
structural integrity and to make the fabrication of the photonic material feasible. 

Such compactness is extremely important if this photonic material shall be used in applications. The omnidrectional 
gap can be exploited, e.g., by using this material as intermediate reflector in tandem solar cells. The problem to be 
solved in such tandem solar cells is to define an intermediate reflector that reflects all light in a certain frequency 
interval and independent of the angle of incidence back into the top cell, leaving the other spectral components 
unaffected. Maximizing the light absorbance in the top cell allows to increase the overall efficiency of the solar cell. 
The proposed photonic material is very promising for such an application, as it acts sufficiently efficient and can be 
fabricated at low costs. 

In conclusion, we have investigated the local density of states for a photonic material made of randomly arranged 
spheres. It was shown that a significant suppression of the LDOS is encountered in frequency domains slightly above 
the Mie resonances of the single sphere. The observation of the effect requires a minimum distance to adjacent spheres. 
The magnitude of the LDOS decays exponentially with a linearly increasing size of the photonic material indicating 
the presence of a photonic band gap of the infinite structure. Covering the spheres with a shell of a permittivity as 
small as feasible allows to observe the same phenomena for self-sustaining ensembles. The theoretical concept might 
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FIG. 4: (color online) (a) Spatial arrangement of a self sustaining ensemble of randomly arranged spheres as generated by the 
sphere dropping method, (b) LDOS for cores with or without shells. The arrangement is only self sustaining if the spheres 
consist of cores and shells. 

be a way to fabricate a photonic material at low costs that suppresses the light propagation along all directions, 
representing an omnidirectional band gap material. 
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